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ABSTRACT 


Shock  wave  attenuation  affects  experimental  results  by  reducing  the 
shock  velocity  as  it  travels  down  the  shock  tube.  Until  recently,  attenuation 
has  not  been  readily  predictable  with  accuracy  and  was  simply  accepted 
during  the  experimental  process.  A  mass  loss  theory  was  developed  by  Dr. 
Harold  Mirels  which  accounted  for  the  shock  wave  attenuation  by 
considering  the  formed  boundary  layer  due  to  the  passage  of  a  shock  wave. 
The  purpose  of  this  thesis  is  to  develop  a  computer  code,  which  would  apply 
Mirels'  theory  to  a  prediction  model,  and  validate  the  results  with  acoustic 
theory  and  Mirels'  data.  It  will  then  be  used  to  determine  an  accurate  shock 
velocity  for  an  additional  shock  tube.  The  computer  results  will  be  proven  to 
match  the  experimental  data  well.  However,  additional  corrections  involving 
friction,  heating,  and  real  gas  effects  are  needed  to  produce  excellent  results. 
It  will  also  be  discovered  that  the  mass  loss  correction  results  are  a 
significant  part  of  the  total  shock  wave  attenuation. 
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CHAPTER  I 


INTRODUCTION 


INTRODUCTION 

Since  its  invention  in  the  early  1900's,  the  shock  tube  has  provided  a 
simple  and  inexpensive  high  speed  flow  phenomena  testing  facility.  Its 
versatility  has  earned  itself  a  regime  of  its  own  in  the  world  of  aerospace 
testing.  Because  of  its  wide  spread  use,  a  large  interest  in  the  shock  tube 
technology  has  ensued,  concerning  the  characteristics  and  performance  of  the 
shock  tube.  These  interests  range  from  attenuation  to  flow  zones. 

Early  on  in  shock  tube  research,  the  attenuation  phenomenon 
interfered  with  accurate  shock  wave  velocity  predictions.  Initial 
investigations  were  accomplished,  but  reliable,  conclusive  results  were  never 
realized.  In  the  1950's  Dr.  Harold  Mirels  arrived  at  a  theory  which  explained 
shock  wave  attenuation  due  to  the  viscous  wall  boundary  layer  interaction 
produced  by  the  passage  of  the  shock  wave.  His  theory  dealt  with  mass 
accountability  in  the  shock  tube  system.  A  mass  transfer  into  the  boundary 
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layer  was  causing  perturbations,  resulting  in  the  attenuation  of  the  shock 
wave  (Mirels,  1956a).  Building  on  the  theory,  he  developed  a  set  of 
equations  which  addressed  the  problem  of  attenuation  in  several  different 
cases. 

Despite  these  developments,  most  laboratories  simply  lived  with  the 
attenuation  problem.  Tedious  calculations  producing  accurate  results  were 
inefficient  when  estimates  provided  enough  information.  Since  shock  tubes 
were  small  and  simple  to  operate  at  the  time,  multiple  runs  to  achieve  a 
single  test  condition  were  not  unmanageable.  The  advent  of  larger  and  more 
complicated  shock  tubes  increased  the  demand  for  better  predictions.  The 
computer  age  allowed  for  better  measurement  techniques.  More  accuracy  in 
measurements  demanded  more  accuracy  in  predictions.  With  this  in  mind, 
the  purpose  of  this  thesis  becomes  apparent. 

PURPOSE 

The  purpose  of  this  thesis  is  to  provided  a  computer  model  of  shock 
tube  flow  characteristics  while  accounting  for  the  shock  wave  attenuation 
based  on  Mirels  theory.  Shock  wave  attenuation  is  made  up  of  several 
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aspects  including:  friction,  heating,  real  gas  effects,  and  Mirels'  mass  loss 
theory.  This  study  will  focus  on  the  mass  loss  due  to  a  turbulent  wall 
boundary  layer  interaction.  The  code  is  intended  to  be  as  versatile  as  the 
shock  tube.  It  should  be  applicable  to  all  shock  tube  configurations  and 
produce  an  accurate  prediction  of  the  shock  wave  velocity,  in  addition  to 
other  important  shock  tube  parameters. 
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CHAPTER  II 


SHOCK  TUBES 


INTRODUCTION 

A  shock  tube  can  be  characterized  as  a  simple  machine  which 
generates  complex  phenomena.  It  consists,  basically,  of  a  rigid  tube  divided 
into  two  sections.  These  sections  are  separated  by  a  gas-tight  diaphragm 
which  is  mounted  normal  to  the  axis.  Initially,  a  significant  pressure 
difference  exists  between  the  two  sections.  The  high  pressure  section  is 
called  the  compression  chamber  while  the  low  pressure  section  is  known  as 
the  expansion  chamber.  At  time  zero,  when  the  diaphragm  is  ruptured,  the 
pressure  begins  to  equalize  in  the  form  of  a  shock  wave  moving  into  the 
expansion  chamber  and  a  rarefaction  wave  moving  into  the  compression 
chamber  (e.g.:  Wright,  1961).  Figure  1  is  a  simple  shock  tube  schematic, 
where  P4  is  the  high  pressure  in  the  expansion  chamber  and  Pj  is  the  low 
pressure  in  the  compression  chamber.  For  the  remainder  of  this  text,  the 
expansion  chamber  will  be  noted  as  the  driver  section,  while  the  compression 
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Pressurization  Valve 


Release  Valve-,^ 


Driver  Section  ^ 

Diaphragm 


Driven  Section 


Figure  1 .  Simple  shock  tube 


chamber  will  be  referred  to  as  the  driven  section.  After  the  diaphragm  has 
been  ruptured,  the  shock  wave  will  move  from  left  to  right,  starting  at  the 
diaphragm.  An  induced  air  flow  will  follow  the  shock  wave.  An  expansion 
wave  will  move  from  right  to  left,  beginning  at  the  diaplu-agm.  There  are 
several  different  versions  and  modifications  which  can  be  made  to  the  simple 
shock  tube.  Only  those  relevant  to  this  text  will  be  mentioned  here. 


HISTORY 

The  shock  tube  appeared  just  prior  to  the  1900's,  when  Vielle  proved 
that  the  shock  wave  moved  faster  than  the  speed  of  sound  down  the  tube. 
After  this,  the  shock  tube  was  dormant  until  the  advent  of  World  War  11. 
Reynolds  used  the  shock  tube  in  1943  to  calibrate  piezoelectric  pressure 
gauges  used  in  blast  wave  measurements.  In  1945,  a  method  for 
photographing  the  position  of  a  shock  front  was  perfected  and  a  detailed 
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study  of  reflected  waves  was  initiated  by  L.  G.  Smith  (Wright,  1961).  After 
the  war,  shock  tubes  seemed  to  become  an  institution  of  their  own.  W. 
Bleakney  developed  a  way  to  determine  pressure  data  for  the  diffraction  of  a 
shock  wave  around  obstacles.  He  used  a  Mach  Zehnder  interferometer  to 
find  the  quantitative  density  which  was  then  used  to  find  the  pressure 
(Wright,  1961).  From  1948  to  1958,  shock  tube  studies  concentrated  on 
one-dimensional  flows  induced  by  shock  and  rarefaction  waves  while  a 
secondary  arena  involved  contact  surface  interactions  (Glass,  1981).  This 
decade  also  saw  the  advent  of  many  problems  with  shock  tube  experiments. 
These  phenomena  hindered  ideal  test  conditions  and  still  present  problems 
today.  The  first  apparent  problem  was  the  diaphragm  burst.  The  ruptured 
diaphragm's  irregular  shape  occupied  some  of  the  tube's  cross-sectional  area. 
In  addition,  barbs  from  the  diaphragm  induced  turbulence  in  flow  through  the 
ruptured  diaphragm,  as  shown  in  Figure  2.  Because  of  this,  flow  quality  was 
degraded  by  the  burst  diaphragm.  The  next  problem  identified  dealt  with  the 
shock  velocity.  It  was  found  that  as  the  shock  moved  down  the  tube,  its 
velocity  and  Mach  number  decreased.  Studies  showed  that  the  rate 
reduction  was  due  to  attenuation  of  the  shock  wave.  Interaction  of  the  side 
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Shock  Tube 
Cross-Section 

Barb 


Figure  2.  Cross-sectional  sketch  of  burst  diaphragm  in  shock  tube 


wall  boundary  layer  induced  by  the  shock  wave  was  responsible  for  the 
attenuation  and  an  acceleration  of  the  contact  surface.  This  phenomena 
affected  characteristics  of  each  run  and  shortened  the  usable  test  time.  An 
additional  problem  caused  by  boundary  layer  interactions  was  flow  gradients 
in  the  driven  section  (Glass,  1981).  The  attenuation  due  to  the  boundary 
layer  interaction  is  the  reason  for  this  thesis. 

Despite  these  shortcomings,  the  shock  tube  had  the  advantage  of  being 
versatile  and  relatively  inexpensive.  Its  appearance  in  university, 
commercial,  and  government  facilities  sparked  the  introduction  of  some 
useful  modifications.  These  included:  shock  tunnels,  multidiaphragm  shock 
tubes,  and  expansion  tunnels  (e.g.:  Glass,  1981).  The  Engineering  Research 
Institute,  at  the  University  of  Michigan  used  a  shock  tube  to  drive  an 
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intermittent  wind  tunnel  (Wright,  1961).  These  hybrids  allowed  the  shock 
tube  to  secure  a  strong  foothold  in  the  fiiture  of  aerospace  studies. 

Work  at  the  University  of  Toronto,  Institute  for  Aerospace  Studies, 
(UTIAS),  focused  on  prominent  characteristics  of  the  shock  tube.  One  of 
these  characteristics  was  wave  interactions.  Shock  wave  refractions  at  a 
contact  surface  and  through  a  gas  layer  were  examined.  It  was  found  that 
special  refracted  shock  waves  would  reflect  only  Mach  wave  from  a  contact 
surface.  The  result  was  "the  reflected  shock  wave,  tailored-interface, 
shock-tunnel  operation,  thereby  providing  a  constant-pressure  high-enthalpy 
reservoir  of  gas  for  expansion  through  supersonic  or  hypersonic  nozzles" 
(Glass,  1981).  Continuing  experiments  dealt  with  protective  gas  barriers  for 
attenuating  blast  waves,  one  shock  wave  overtaking  another,  and  speed  of 
sound  measurements  for  different  gases.  Shock  wave/shock  wave  and  shock 
wave/rarefaction  wave  collisions  were  also  investigated.  It  was  found  that 
very  high  temperatures  could  be  reached  through  the  collision  of  two  shock 
waves  of  equal  strength.  The  last  important  point  in  this  wave  of 
experiments  at  UTIAS  was  the  finding  that  strong  shock  wave  predictions 
had  to  include  real  gas  effects  in  their  calculations  in  order  to  match  the 
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experimental  results.  This  was  probably  no  surprise,  since  strong  shock 
waves  will  induce  high  enthalpy  flows.  However,  it  did  imply  that 
one-dimensional  shock  theory  was,  at  best,  and  estimate  in  most  experiments 
involving  shock  tubes. 

In  the  1960's,  the  shock  tube  applied  itself  outside  the  world  of 
aerodynamics  by  producing  industrial  grade  diamonds  from  graphite.  The 
pressures  generated  to  form  the  diamonds  were  measured  in  the  megabar 
range,  while  temperatures  soared  above  the  million  degree  mark.  Although 
diamond  sizes  were  small,  (10-20  micrometers),  several  companies, 
including  3M  Canada  Ltd.,  continued  the  study  with  large  success.  This 
configuration  also  developed  fusion  plasmas  in  deuterium  (Glass,  1981). 
Another  use  of  the  shock  tube  outside  the  aerospace  industry  is  the 
development  of  high  power  gasdynamic  and  chemical  lasers  (Anderson, 
1990). 

An  example  of  a  hypersonic  conditions  simulator  is  the  G-Range 
Impulse  Facility  at  the  Arnold  Engineering  and  Development  Center 
(AEDC).  This  facility  doubles  as  a  high  performance  free  piston  shock 
tunnel  and  a  two-stage  light  gas  gun.  Predictions  for  the  shock  tunnel  include 
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stagnation  pressures  of  10,000  atm.  and  enthalpies  of  20  MJ/kg  (Maus, 
Laster,  and  Homung,  1992).  Figure  3  shows  the  launcher  and  shock  tunnel 
configurations.  The  shock  tunnel  offers  the  advantage  of  a  recoverable 
model  while  still  maintaining  high  Mach  numbers  and  temperatures.  These 
facilities  represent  the  latest  in  hypersonic  testing  mechanisms. 

The  shock  tube  has  become,  and  will  probably  remain  a  mainstay  of 
the  aerospace  testing  community.  It  is  an  inexpensive  resource  which  can  be 
found  throughout  the  world  in  government,  commercial,  and  university 
laboratories.  Because  of  its  adaptability  and  versatility,  it  provides  a  variety 
of  test  conditions,  including  those  conditions  which  cannot  be  duplicated  by 
any  other  type  of  facility.  Derivatives  of  the  shock  tube  technology  have 
provided  results  in  areas  which  were  previously  unsolved.  The  shock  tube 
has  even  expanded  its  usefulness  into  other  industries.  Despite  reductions  in 
funding,  advances  in  shock  tubes  and  testing  techniques  continue  to  set  the 
stage  for  more  practical  results  and  widespread  applications. 
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COMPRESSION 
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SECOND  DIAPHRAGM 
inCATPD  AT  THROAT 

Figure  3.  Impulse  Facility  Launcher  and  shock  tunnel  at  ALUC 
(Mans,  Faster,  and  Homung) 


SHOCK  TUBE  THEORY 


The  basis  for  shock  tube  theory  is  the  assumptions  and  derivations  for 
a  moving  normal  shock  wave.  In  certain  cases,  including  this  thesis,  the 
expansion  wave  is  irrelevant  because  testing  is  completed  by  the  time  the 
expansion  wave  has  interfered.  Therefore,  the  expansion  wave  process  is 
excluded  from  the  derivations  of  moving  normal  shock  wave  equations  for 
shock  tubes.  Derivation  of  these  equations  can  be  found  in  Anderson's 
Modem  Compressible  Flow,  and  begins  with  the  following  forms  of  the 
continuity,  momentum,  and  energy  equations  for  a  shock  fixed  coordinate 
system: 

PlUi  =  P2U2  (2.1) 

Pi  +  P  1^1^  =  P2  +  9  2^2  (2-2) 

hj  +  UjVl  =  h2  +  U2V2  (2.3) 

where  Uj  is  the  velocity  of  the  gas  ahead  of  the  shock  wave  and  U2  is  the 

velocity  of  the  gas  behind  the  shock  wave.  Both  velocities  are  relative  to  the 

shock  wave  (Anderson,  1990).  Substituting  W  for  u,  and  W-Up  for  U2  yields 

the  following  set  of  equations; 

^2  -  61  =  (Vi-V2)(P,  +  P2y2  (2.4) 
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Equation  2.4  is  the  Hugoniot  equation,  which  is  the  same  for  moving  and 
stationary  shocks.  Here  v  =  RT/p  and  e  =  c^T.  The  rest  of  the  relationships 


follow  for  the  moving  shock  wave  in  laboratory  coordinates: 

M,  =  W/ai 


£2 

Pi 


f  Zli+El  ^ 

Y-l  Pi 

V  1+^^  J 

V  Y-l  Pi  y 


I  Y+l  P2 
£2  _  Y-l  Pi 
Pi  ~  Y+1  P2 

Y-l  Pi 


(2.5) 


(2.6) 


(2.7) 


U 


P 


_  £1  I 

”  T  VPiy 


Y+l 


Y-l 

Y+l 


(2.8) 


(2.9) 

where  the  specific  heat  ratio  is;  y  =  CJC^.  Figure  4  is  reproduced  from 
Anderson  and  shows  the  one-dimensional  theoretical  pressure  distribution 

throughout  a  shock  tube.  The  shock  Mach  number,  M^,  is  found  by  dividing 

W  by  a^.  The  shock  tube  is  labeled  with  notation  corresponding  to  the  above 

equations.  The  shock  wave  and  contact  surface  move  from  left  to  right  while 
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x-Distance 


Pi 

Figure  4.  Shock  tube  flow  with  theoretical  pressure  distribution  after 
diaphragm  burst  (Anderson,  1990) 

the  expansion  wave  moves  from  right  to  left.  Regions  1  and  4  are 
intentionally  identified  because  they  will  remain  in  their  initial  state 
throughout  the  test.  Region  2  consists  of  the  driven  gas  after  shock  wave 
interaction.  Region  3  is  the  most  complicated  region  since  it  includes  both 
the  driver  and  driven  gases,  should  they  be  different,  and  has  interacted  with 
the  shock  wave  and  the  contact  surface.  For  most  testing  and  computations 
however,  this  region  is  ignored  due  to  its  small  role  in  the  shock  wave 
propagation. 

An  interesting  note  concerning  the  flow  velocity  behind  the  shock 
wave  is  its  limiting  value.  If  a  calorically  perfect  gas  is  assumed  with  a  y  of 
1 .4  and  a  pressure  ratio  approaching  infinity,  the  maximum  Mach  number  of 
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the  flow  behind  the  shock  wave  is  1.89  (Anderson,  1990).  This  can  be 
confirmed  by  taking  the  limit  of  the  equation  for  the  ratio  of  Up/a2  as  shown 
in  equation  7.19  of  Anderson.  This  unsteady  flow  theory  forms  the  basis  for 
shock  tube  analysis  and  calculations. 
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CHAPTER  III 


MIRELS’  SHOCK  WAVE  ATTENUATION  THEORY 

INTRODUCTION 

One-dimensional  shock  tube  theory  gives  a  good  estimate  for  the 
shock  Mach  number  in  a  shock  tube.  For  low  Mach  numbers  (M  <  1.5 
approximately),  the  Mach  number  estimate  may  even  be  considered  accurate. 
According  to  classical  theory,  this  shock  Mach  number  is  considered  a 
constant  from  the  diaphragm  to  the  end  of  the  driven  section  for  a  constant 
area  tube.  As  stated  in  Chapter  2,  it  was  found  that  the  shock  wave  actually 
attenuated  with  distance,  decreasing  the  shock  Mach  number  significantly. 
This  problem  led  to  studies  in  the  1950's  by  Dr.  Harold  Mirels,  who 
concentrated  his  efforts  on  a  more  accurate  prediction  of  the  available  shock 
tube  test  time.  His  theory  is  based  on  the  development  of  the  boundary  layer 
in  the  shock  tube  (Mirels,  1956a).  The  shock  tube  boundary  layer  was  never 
included  in  elementary  theory  due  to  the  inviscid  flow  assumption. 
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BASIC  THEORY 


The  introduction  of  the  boundary  layer  into  the  theoretical  predictions 
presents  some  revelations.  To  proceed  with  his  analysis,  Mirels  assumes  that 
the  boundary  layer  is  thin  compared  with  the  shock  tube  diameter.  The 
assumption  is  based  on  the  premise  that  most  shock  tubes  are  designed  for  a 
uniform  core  of  potential  flow.  The  thin  boundary  layer  forms  along  the 
walls  of  the  shock  tube  behind  the  shock  wave  as  shown  in  Figure  5. 
According  to  Mirels'  theory  this  boundary  layer  emits  weak  expansion 
waves.  These  acoustic  waves  over  take  the  shock  wave,  causing  its 
attenuation  (Mirels,  1956a).  Previous  theories  for  the  shock  wave 
attenuation  are  addressed  by  Mirels,  but  will  not  be  discussed  here. 


Expansion  Wave  Shock  Wave 


Figure  5.  Shock  tube  flow  with  boundary  layers  (Mirels,  1956) 
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Mirels'  analysis  (1956a)  begins  with  the  assumption  that  the  motion  in 
the  shock  tube  can  be  considered  one-dimensional  unsteady  flow  through  a 
constant  area.  Each  incremental  cross-section  of  the  tube  is  considered  as  a 
weak  mass  source  which  perturbs  flow  through  the  tube  through  the 
generation  of  pressure  waves.  Mirels  begins  his  derivations  with  three  basic 
equations.  These  equations  are; 


Momentum; 

(  5Aw  ,  dAw  1 

PV^  +  “-arj  = 

dAp 

dx 

(3.1) 

Continuity: 

SAp  ^  ^BAu  ,  dAp 

at  +  p  &  +  “  at 

=  m 

(3.2) 

Isentropy; 

Ap  =  fl^Ap 

(3.3) 

Here  m  =  m(x,0,  which  constitutes  the  rate  of  mass  addition  per  unit 
cross-sectional  area  per  unit  x.  The  equations  for  the  net  perturbation  at  any 
point  X,  t  are  expressed  as  follows: 

The  +  superscript  corresponds  to  a  positive  x-direction  wave  motion  and  the 
-  superscript  associates  with  a  negative  x-direction  motion.  The  variable  ^  is 
an  integration  variable  for  x. 
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These  perturbation  equations  are  fundamental  to  acoustic  wave  theory. 
They  will  be  recalled  in  the  following  chapter  to  assist  in  basic  computer 

code  validation.  The  next  step  Mirels  takes  is  to  relate  the  ratio  of  the 

change  in  pressure  over  the  total  pressure  to  the  normal  velocity,  v.  A 
positive  normal  velocity  yields  a  compression  wave  while  a  negative  velocity 

forms  an  expansion  wave,  where  the  term  "positive"  in  this  case,  applies  to 

an  inward  normal  velocity  away  from  the  wall. 

For  the  derivation  of  the  general  shock  attenuation  formulas,  flow  in 

the  shock  tube  is  assumed  to  include  ideal  basic  flow  and  small  perturbations 

from  the  boundary  layer.  The  ideal  flow  expansion  wave  is  assumed  to  have 

a  thickness  of  zero  with  a  speed  equivalent  to  the  speed  of  sound  in  region  4, 

shown  in  Figure  5.  These  derivations  are  not  relevant  here,  but  can  be  found 

in  NACA  Report  1333. 

At  this  point,  Mirels'  theory  departs  in  different  directions.  One 
direction  addresses  wholly  laminar  flow  in  a  specified  region,  while  another 
solves  the  wholly  turbulent  flow  region.  An  obvious  question  concerns  the 
case  where  a  transition  from  laminar  to  turbulent  occurs  within  a  single 
region  of  the  shock  tube.  Although  Mirels  does  provide  a  method  to  address 
this  case,  it  is  difficult  to  apply.  In  addition,  acceptable  methods  for 
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determining  the  flow  transition  point  have  not  been  developed  (Mirels, 
1956a).  Therefore,  this  case  will  not  be  treated  here. 

For  the  standard  shock  tube  with  a  driver  section,  represented  as 
region  3  in  Figure  5,  and  a  driven  section,  shown  as  region  2,  four  different 
cases  are  apparent.  The  first  is  the  wholly  laminar  case,  where  both  region  2 
and  region  3  contain  laminar  flow.  The  second  and  third  cases  concern  a  mix 
of  flow,  where  one  region  is  wholly  turbulent  and  one  is  wholly  laminar.  The 
last  occurs  when  both  region  2  and  region  3  are  wholly  turbulent.  This  last 
situation  is  the  most  common  and  the  only  case  addressed  in  this  thesis.  The 
reasoning  behind  this  decision  is  simple.  The  developed  code  from  this 
report  will  first  be  applied  to  shock  tubes  located  at  Arnold  Engineering  and 
Development  Center,  specifically,  the  G-Range  high  performance  free  piston 
shock  tunnel.  The  majority  of  these  shock  tubes  operate  in  the  high 
Reynolds  number  range  and  therefore  exclusively  involve  wholly  turbulent 
flow  in  both  the  driven  and  driver  sections  of  the  shock  tubes. 

As  mentioned  before,  normal  velocities  occur  as  a  result  of  the  mass 
loss  effects.  Other  effects  noted  by  Mirels  concern  dissipation  and  heat 
transfer.  Dissipation  increases  the  normal  velocity.  Heat  transfer  increases 
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the  normal  velocity  when  passing  from  wall  to  boundary  layer.  The  opposite 
situation  has  a  decreasing  effect  on  the  normal  velocity. 

TURBULENT  BOUNDARY  LAYER 

Integral  methods  are  used  to  solve  the  turbulent  boundary  layer 
problem.  This  solution  will  be  addressed  in  detail.  Mirels  (1956b)  used 
empirical,  semi-infinite  flat-plate,  boundary  layer  data  applied  to  a  moving 
wall  to  develop  these  solutions.  The  ensuing  set  of  derivations  follows  the 
development  in  NACA  TN  3712.  For  a  fixed  shock  coordinate  system,  the 
steady  flow  fluid  dynamic  equations  apply  in  the  region  bounded  by  the 
shock  and  the  interface.  Since  one-dimensional  shock  tube  theory  predicts  a 
constant  pressure  in  this  region,  the  boundary  layer  equations  can  be  written 
as; 

d^u  dpv  ^ 

Continuity:  ^ 

Momentum;  +  pv^  =  ^  (3.7) 

Energy: 

with  the  boundary  conditions: 
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u{x,  0)  =  Um>  u{x,  Co)  =  Ue  v(jr,  0)  =  0 
it  is  possible  to  find  the  integral  form  of  the  momentum  equation.  Figure  6 

on  the  following  page  shows  the  fixed  shock  coordinate  system  used  for  the 

derivation.  Here,  W  fi-om  Equation  2.9  is  equal  to  U^,  and  u^  =  -  U2. 

Integrating  the  momentum  equation  with  respect  to  y  fi*om  zero  to  5  yields: 


io  P’^%dy  +  \l(>'^%dy  =  -T 


Rearranging  the  continuity  equation  provides: 


Substitution  of  equation  3.10  into  3.9  gives: 


(3.9) 

(3.10) 


(3.11) 


Integration  by  parts  allows  for  the  following: 

-Ue\l^-^<fy  +  \l(pU^  +  U^)dy  =  -T„  (3.12) 

Combination  of  like  terms  and  simplification  confirm  the  integral  form  of  the 


momentum  equation  as: 

f5  pM  _  jLx  J 
n  ~  Jo  PeMef*  ~  dx 

Pe^e 


(3.13) 


The  term  u/u^  is  the  average  velocity.  In  order  to  continue,  some 


assumptions  and  definitions  must  be  noted.  First,  the  boundary  layer 
thickness  will  be  noted  as  8  and  a  similarity  parameter  =yl%.  To  express 
u/Ug  as  a  function  of  the  similarity  parameter,  Mirels  assumes  the  average 
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Shock  uw 

Figure  6.  Fixed  shock  coordinate  system 


velocity  will  behave  according  to  the  seventh-power  law,  which  was  derived 
for  a  turbulent  boundary  layer  on  a  semi-infinite  flat  plate  according  to 
Mirels.  Therefore, 

^  =  (^7  for  0<(^<1  (3.14) 

This  format  has  provided  good  results  for  the  entire  range  of  Mach  number 

flows.  The  seventh-power  law  also  works  well  for  skin  friction  data  (Mirels, 

1956b).  To  include  the  moving  wall  profile,  velocities  relative  to  the  wall 


necessary  calculations  include  the  following  ratios  with  appropriate 


coefficients  (Mirels,  1964): 
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^  =  i  =  +  for  0<C^1  (3.17) 

^  =  ^=1  for  1<(;  (3.18) 

Where  b  =  (h/h^)-!  and  c  =  h/h^-hg  o/h^..  Mirels  combines  the  previous  two 
sets  of  equations.  Equations  3.15  to  3.18,  to  find  the  boundary  layer 


momentum  thickness; 


0 


6  ■  T’m- 

and  the  boundary  layer  displacement  thickness: 

5*  1  nTe 


= 7^(1  -  + (1  -  2^)h  -  (1  -  s)/8]  (3.19) 


=  l-7^[^/6  +  (l-t)/7] 


The  L,  terms  are  defined  as; 


1  z^dz 
0  \+bz-cz^  ' 


(3.20) 


(3.21) 


The  next  step  Mirels  takes  is  to  work  toward  boundary  layer  relationships. 
Assuming  the  wall  temperature  is  constant,  0/6  is  independent  of  x.  The 


integral  form  of  the  momentum  equation  3.13,  can  be  rewritten. 


Tw  _ 

PeMe  5  dx 

Here,  the  Blasius  relation  is  required  to  associate  Xw 


(3.22) 


and  5  which,  for 


incompressible  turbulent  flow  on  a  semi-infinite  flat  plate  is: 


^  =  0.0225(f^)  j 

pgMe  WeO 


(3.23) 
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To  solve  for  the  boundary  layer  thickness,  a  new  subscript  "m"  is  introduced 
in  the  following  equations  and  denotes  a  mean  value  for  the  appropriate 


variable. 


=  0.0225(-^^)J 

QmUe  "eO 


(3.24) 


^  =  0.0225(K;f4)J 

QeUe  ^eO 


(3.25) 


rf  A  r 

)  \Tr„) 


(3.26) 


An  estimate  for  the  mean  temperature  value  can  be  found  with  the  equation 
=  0.5(T^+TJ+0.22(Tr-TJ  (Mirels,  1956b).  Multiplying  a  simple  factor  to 
the  Equation  3.25  above  gives; 


9e 


7-^4—  =  0.0225(|)(-^^)Jf^ 

(Ue-Uw)^  ^^\Ue-Uyv\^‘'  \Ue-llw\ 


(3.27) 


Simplifying  this  equation  yields; 

^  =  0.0225^1 (3.28) 

Another  substitution  involving  Equation  3.22  provides  the  following; 

3  1 


If  =  0.0225(^(1 -^)|  1-^1 


(3.29) 


If  is  considered  constant,  the  Equation  3.29  can  be  integrated,  giving; 


5  =  0.0574((t)-!-^)?|l 


0/6 


Ue 


(3.30) 
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In  addition  to  these  equations,  some  relationships  dealing  with  enthalpy  are 


required.  The  mean  static  enthalpy  is  defined  by: 

^  =  0.5(^+l)  +  0.22(^-l)  (3.31) 

where  h^  is  the  recovery  enthalpy  (Mirels,  1964).  The  previous  ratios  are 


defined  for  an  ideal  gas  as: 


hr 

hefi 


1  + 


(W-lf 

ZW-l 


and 


hw 

he,0 


W(Z-W) 

ZW-l 


(3.32,  3.33) 


where 

and  and  KO)  =  (Pr)’ 

The  r(0)  term  is  the  recovery  factor.  For  a  Y=1 .4,  r(0)  is  found  to  be  0.897 
and  for  '7=1.667,  r(0)  =  0.875  (Mirels,  1964). 

Although  Mirels  provides  additional  equations  for  use  in  many 


situations,  only  those  listed  here  are  used  in  the  prediction  code  to  reproduce 


Mirels'  approach  to  shock  wave  attenuation.  Integrals  were  performed  using 


Simpson's  Rule.  All  other  equations  were  combined  in  a  single  subroutine 


which  was  then  attached  to  an  existing  computer  code. 
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CHAPTER  IV 


SHOCK  TUBE  PREDICTION  PROGRAM 

INTRODUCTION 

A  working  computer  program  was  the  quantifying  objective  of  this 
thesis.  The  code  would  be  modified  for  each  individual  shock  tube.  Inputs 
would  be  allowed  for  various  runs.  Calculations  would  include  shock  tube 
characteristics  and  other  relevant  data  for  the  test.  The  difference  between 
this  prediction  code  and  those  already  in  place  would  be  the  accurate 
determination  of  the  shock  Mach  number,  shock  tube  characteristics, 
reflected  shock  values,  and  other  appropriate  variables  formerly  estimated. 

Although  optimistic,  expected  results  were  not  anticipated  to  be  1 00% 
accurate.  Instead,  the  mass  loss  attenuation  theory  application  was 
considered  a  significant  step  in  the  right  direction  to  solving  the  continuing 
problem  of  the  over-predicted  shock  velocity.  With  a  solid  foundation  for 
more  accurate  prediction  methods,  refinements  could  be  pursued,  enhancing 
the  effectiveness  of  the  predictions.  This  would  in  turn  allow  more  efficient 
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testing  in  the  form  of  increased  test  precision,  fewer  required  runs,  and  better 
computational  methods.  The  indirect  advantage  is  monetary  savings,  which 
is  a  welcomed  result  in  any  field  of  endeavor. 

THE  PROGRAM  DESCRIPTION 

The  basic  code  used  for  this  Fortran  computer  program  was  adapted 
fi-om  an  existing  prediction  code  for  the  free  piston  shock  tunnel  located  at 
AEDC.  The  original  code  was  applied  to  the  fi'ee  piston  shock  tunnel  in 
1989  (Mans,  1994).  Despite  its  origins,  the  computer  algorithm  can  be 
applied  to  any  turbulent  boundary  layer  shock  tube. 

This  program  operates  in  increments  of  position  within  increments  of 
time  in  a  Lagrangian  format.  A  tube  is  separated  into  the  two  basic  shock 
tube  sections:  the  driver  section  and  the  driven  section.  Each  has  its  own  set 
of  initial  conditions  corresponding  to  the  appropriate  shock  tube  dimensions 
and  setup.  Each  shock  tube  section  is  divided  into  a  specified  number  of 
zones  which  remains  fixed  for  a  given  run.  These  zones  contain  a  initial 
amount  of  mass,  determined  by  the  volume  of  the  tube  and  type  of  gas  in  the 
tube.  The  mass  is  considered  to  be  concentrated  in  two  halves  at  the  center 
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of  each  respective  volumetric  dividing  line.  Figure  7  depicts  some  of  the 
notation  and  concepts  of  the  code  methodology.  Newtonian  d)mamics 
govern  the  movement  of  the  mass  particles  while  the  basic  theory  solves 
variables  such  as  pressure,  volume,  and  energy.  For  each  cycle,  increments 
of  time  and  distance  are  added  to  advance  the  shock  tube  system.  The 
program  is  set  up  to  terminate  at  the  end  of  a  selected  number  of  cycles  or 
when  the  shock  reaches  the  end  of  the  tube,  whichever  occurs  first.  The 
program  can  be  modified  to  terminate  for  other  events  as  well  or  continue 
into  the  reflected  shock  regime,  although  attenuation  calculations  would 
require  some  modifications. 

Within  a  time  step,  the  program  is  governed  by  five  basic  equations. 
The  first  three  involve  Newton's  second  law.  Velocity  and  displacement  are 
determined  for  each  mass  point  by  these  three  equations; 


^-§  =  iPQj_L-PQj,L-DPF)^ 

(4.1) 

(4.2) 

=  X^j  +  dt^ 

(4.3) 

The  term  PQ  is  a  pressure  term  where  P  is  the  pressure  and  Q  is  a  dissipation 
factor.  Q  is  proportional  to  the  square  of  the  local  velocity  gradient  whose 
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Figure  7.  Mass  points,  cells  and  regions  (Maus,  1991) 


value  is  only  significant  near  the  shock  wave.  DPF  represents  a  drop  in 
pressure  due  to  fiiction  forces  on  the  mass  point  (Mans,  1991). 

The  next  two  equations,  energy  and  equation  of  state,  solve  the  change 
in  energy  and  the  pressure.  The  equation  forms  used  in  the  program  are: 


dE  =  dQw  -  PQdV 

(4.4) 

P=j{E,V)=^^ 

(4.5) 

The  term  dQ^  is  found  using  Newton's  law  of  cooling.  With  these  two 

variables  known,  the  rest  of  the  thermodynamic  variables  can  be  determined. 

The  mass  loss  subroutine,  labeled  SINK,  contains  the  equations  for  the 

attenuation  algorithm  which  were  noted  in  Chapter  3.  The  three  necessary 

integral  terms  are  solved  using  Simpson's  Rule.  The  equations  for  the 

0  5* 

boundary  layer  momentum  thickness,  ~  (Equation  3.19),  and  (Equation 
3.20)  the  displacement  thickness,  help  provide  the  following  relationships. 

5  cc  ^  (4.6) 

This  equation  is  derived  fi'om  Equation  3.30.  The  normal  velocity  is  found  to 


be. 


y=Ue 


db* 

dx 


The  mass  flow  into  the  boundary  layer  is, 

m  =  PeUe^P^X 


(4.7) 


(4.8) 
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This  mass  flow  term  represents  the  net  change  in  mass  leaving  the  potential 
flow  and  entering  the  boundary  layer,  causing  the  perturbations.  Since  this 
mass  flow  term  comes  from  one  unit  of  mass  it  must  be  split  in  half.  Recall 
that  each  mass  point  is  separated  into  two  equal  parts,  residing  in  adjacent 
volume  cells.  Two  different  half  mass  terms  are  used  to  account  for  the  total 
mass  in  the  cell.  To  summarize,  a  mass  flow  into  the  boundary  layer  was 
calculated.  The  mass  traveling  into  the  boundary  layer  was  subtracted  from 
the  mass  in  the  potential  flow,  giving  a  new  half  mass  term.  This  new  half 
mass  term  is  carried  through  on  the  next  iteration  and  used  to  solve  the  next 
increment  of  shock  tube  characteristics. 

METHOD  VALIDATION 

In  order  to  determine  whether  the  computer  code  algorithm  was 
accurate,  acoustic  theory  was  applied  and  solutions  were  compared.  The 
code  was  simplified  to  act  on  a  single  region  tube.  The  initial  state  variables 
were  constant  throughout  the  tube.  At  some  time,  t,  the  mass  of  a  cell  in  the 
center  of  the  tube  was  altered  and  the  resulting  motion  was  examined  and 
compared  with  an  exact  solution  from  acoustic  wave  theory. 
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The  computer  code  was  modified  by  the  author  to  produce  the 
acoustic  results.  The  following  calculations  and  figures  were  provided  for 
comparison  (Mans,  1994).  The  equation  for  the  propagation  of  an  acoustic 
wave  is: 

Ap+(x,0  =  §!!<„ (4.1) 

where,  m{x,t)  =  M^,  for  jc,  <  jc  <  jCj  and  ^  >  0.  Equation  4. 1  is  another  form  of 
Equation  3.4  from  Mirels'  report  (Mirels,  1956a).  The  m  term  is  the  rate  of 
mass  addition  per  unit  cross  sectional  area  per  unit  length.  Figure  8  displays 
the  mass  source  region  and  a  random  point,  x,  for  an  acoustic  wave 
propagation  tube.  Solutions  for  given  intervals  of  t  and  a  fixed  field  point,  x, 
are  as  follows: 


o<?<^ 

0  =  0 

X-X2  ,  X-X 1 

a  <t<  a 

t)  =  ^Ma(x2  -{x-  at)) 

t  > 

1  ^  a 

Ap^{x,  t)  =  fMafc  -x\) 

Source  Region 


X,  X2 


X 

Field  Point 


Figure  8.  Acoustic  wave  tube  diagram  (Maus,  1994) 
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Figure  9.  Fixed  field  point  graph  (Maus,  1994) 

Figure  9  is  a  graphical  summary  of  the  set  of  solutions  for  the  fixed  field 

point.  Here,  =  (jc  -  and  ^2  =  (-^  ■  x.^la.  The  next  set  of  solutions  are 

given  over  the  interval  x  for  a  fixed  time,  t.  Figure  10  is  the  graphical 

representation  for  the  fixed  time  case. 

X>X2  +  Clt  Ap  =  0 

x\+ at  <x  <X2  + at  =  ^Ma{x2-{x- at)) 

x<x\+at  l^p  =  ^Ma{X2-X\) 

With  this  ground  work,  numbers  were  substituted  for  variables  in  the  exact 

solution  and  the  computer  code.  The  following  variables  were  set  at:  length 
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=  1000  cm,  AXq  =  5  cm,  X2  -  jCj  =  25  cm,  diam  =  8  cm,  A  =  50.26  cm^,  F  = 
0.01  (reduction  in  mass  per  step),  a  =  347  m/sec,  and  is  found  to  equal 
0.07958  gm/sec/cm^  These  numbers  yield:  Ap  =  0.03451  bar  and  p  = 
1.0345  bar. 

Figures  11  and  12  show  the  program  generated  results  along  with  the 
exact  solutions  at  t  =  2  and  4  msec,  respectively.  Figure  11  shows  no 
difference  between  the  two  methods.  Figure  12  indicates  a  slight  difference. 
This  difference  increases  with  time,  but  remains  within  the  bounds  of 
acceptable  error.  Therefore,  the  conclusion  can  be  drawn  that  the  program 
models  acoustic  propagations  accurately  indicating  that  it  should  model  the 
mass  loss  algorithm  adequately. 

The  computer  code  has  been  shown  to  match  the  acoustic  theory  well. 
However,  this  evidence  is  not  conclusive.  In  order  to  examine  the  program 
effectiveness,  a  plot  of  the  pressure  distribution  is  necessary.  The  plot 
should  resemble  the  theoretical  plot  in  Figure  4.  Figure  13  shows  an 
acceptable  match  between  the  theory  and  the  computer  results.  Classical 
theory  results  are  also  shown  in  Figure  13. 
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Figure  1 1 .  Acoustic  wave  propagation  for  t  =  2  ms  (Mans,  1994) 
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Acoustic  wave  propagation  for  t  =  4  ms  (Mans,  1994) 


Classical  Pressure  Distribution,  t  =  2.59ms 
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Figure  13.  Pressure  distribution  from  computer  code  for  P4/P1=4.061 


CHAPTER  V 


RESULTS 


INTRODUCTION 

Results  from  the  computer  code  were  easy  to  produce.  Difficulty 
came  when  the  search  was  begun  for  reliable  shock  tube  data  which  included 
the  shock  tube  initial  settings,  dimensions,  and  a  measured  characteristics  at 
known  x  locations.  Additionally,  comparisons  were  made  within  a  single 
run,  comparing  the  mass  loss  algorithm  with  the  classical  approach  and 
examining  the  attenuation  through  graphs. 

Two  shock  tube  configurations  were  selected.  The  first  set  of  cases 
addresses  the  four  runs  used  by  Mirels  in  his  validation  (Mirels,  1956a). 
Here,  the  shock  Mach  number  range  is  approximately  1 .3  to  1.8.  The  second 
set  involves  data  from  a  research  shock  tube  for  diagnostic  development  at 
Arnold  Engineering  and  Development  Center.  The  Mach  number  for  this  set 
of  data  is  approximately  4.0.  The  results  varied  from  excellent  to  fair  within 
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each  case.  Adjustments  to  the  code  were  made  as  a  consequence  of  initial 
figures. 

RESULTS 

The  attenuated/ideal  pressure  ratio,  P2,d^2’  has  been  plotted  by  Mirels 
for  four  separate  shock  tube  runs.  The  ratio  represents  an  excellent  indicator 
of  the  shock  wave  attenuation  in  a  shock  tube.  For  this  reason,  the  data  was 
used  to  compare  with  results  from  the  computer  code.  The  code  results 
were  compared  with  Mirels'  graphs  to  ensure  conformity.  The  computer 
code  graphs  are  included  here  for  each  of  the  four  pressure  ratios.  In  the 
following  figures,  as  in  Mirels'  graphs,  ?2j  represents  the  calculated  pressure 
immediately  behind  the  shock  wave  using  the  mass  loss  equations.  In  order 
to  maintain  uniformity,  ten  pressure  terms  behind  the  shock  were  averaged  to 
find  P2d-  ^2  the  pressure  behind  the  shock,  calculated  by  the  computer 
code  without  the  mass  loss  corrections.  It  was  also  determined  by  the 
average  of  ten  pressure  terms  immediately  behind  the  shock  wave.  Figures 
14  through  17  show  the  results  for  the  four  initial  pressure  ratios.  The 
squares  indicate  the  experimental  measured  points  from  Mirels'  source  data. 
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Figure  14.  Comparison  graph  generated  from  computer  program  for  M=1.313 
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Figure  15.  Comparison  graph  generated  from  computer  program  for  M=1 .406 
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Comparison  graph  generated  from  computer  program  fr 
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Figure  17.  Comparison  graph  generated  from  computer  program  for  M=1 .718 


P  ratio  is  the  calculated  computer  results.  The  line  through  the  P  ratio  points 
is  a  polynomial  regression.  The  line  containing  the  triangular  flags  shows 
Mirels'  calculated  results  (Mirels,  1956a).  The  Mach  numbers  are  averaged 
over  the  length  of  the  tube  and  are  calculated  with  the  mass  loss  subroutine. 
Although  the  computer  program  solutions  are  slightly  erratic,  they  are 
acceptable.  The  trend  for  all  four  graphs  is  an  overestimated  pressure  ratio, 
indicating  an  underestimated  attenuation.  However,  a  difference  between  the 
results  of  the  mass  loss  algorithm  and  the  classical  solutions  is  evident. 
Uncorrected  Mach  numbers  calculated  by  the  computer  program  for  the  four 
previous  cases  were;  M  =  1.34,  1.45,  1.52,  and  1.78,  respectively. 

To  demonstrate  the  attenuation  phenomenon,  plots  were  generated 
comparing  the  shock  wave  locations,  with  respect  to  time,  for  program  runs 
with  and  without  the  mass  loss  subroutine.  Figure  18  is  one  of  those  plots, 
showing  the  mass  loss  data  as  diamonds  (upper  line)  and  the  classical  data  as 
a  vertical  dash  (lower  line).  This  data  is  based  on  inputs  and  characteristics 
from  the  diagnostic  development  shock  tube  at  AEDC.  This  shock  tube  uses 
a  Helium  driver  and  a  Nitrogen  driven  section.  For  the  case  in  Figure  1 8,  the 
driver  pressure,  P4,  is  860  psi.  The  driven  section  pressure,  P,,  is  1.47  psi. 
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Figure  18.  Shock  wave  attenuation  comparison  graph 


The  mass  loss  line  in  Figure  18  obviously  slows  with  time  and  distance. 
Plots  for  other  cases,  including  Mirels',  show  similar  trends.  The  undulations 
in  the  plots  illustrate  the  difficulty  in  determining  the  location  of  a 
computationally  captured  shock  wave. 

The  next  set  of  graphs  deal  with  the  diagnostic  development  shock 
tube  at  AEDC.  Three  cases  are  examined  and  labeled:  A,  B,  and  C.  All 
three  cases  have  a  driver  pressure  (PJ  of  860  psi.  The  driven  pressures  (Pj) 
for  the  runs  are:  76  torr  (1.47  psi),  50  torr  (0.967  psi),  and  115  torr  (2.22 
psi),  respectively.  Each  case  was  computed  with  and  without  the  mass  loss 
subroutine.  Both  shock  speed  calculations  are  plotted  in  the  following 
figures,  as  are  the  measured  shock  speed  data  (Smith,  1994). 

Figure  19  shows  Case  A.  Excellent  agreement  is  observed  between 
the  mass  loss  prediction  values  and  the  measured  values.  Computer  results 
without  the  mass  loss  subroutine  differ  by  almost  30%.  Unfortunately,  this 
case  cannot  stand  alone  and  Figure  20  displays  over-predicted  results  for 
Case  B.  Case  B  contains  the  largest  initial  pressure  ratio  of  the  tliree  runs.  It 
could  therefore  be  assumed  that  Case  C,  having  the  lowest  pressure  ratio, 
will  have  under-predicted  results.  Figure  21  proves  this  conclusion  is 
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Figure  19.  Diagnostic  development  shock  tube  shock  wave  velocity:  Case  A 


Case  B 


Diagnostic  development  shock  tube  shock  wave 


Case  C 
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Diagnostic  development  shock  tube  shock  wave 


incorrect  although  the  initial  tendency  favored  the  hypothesis. 

One  last  comment  concerns  a  comparison  between  mass  loss 
corrections  and  heat/friction  corrects.  The  first  case  fi-om  Mirels,  P4/P,  = 
4.061,  was  run  in  four  different  configurations  (Mirels,  1956a).  The  first  run 
included  mass  loss  and  heating  and  fiiction  corrections.  The  average  shock 
Mach  number  for  this  run  was  calculated  at  1.31.  When  only  mass  loss 
corrections  were  considered,  the  new  Mach  number  was  1.33.  Conversely, 
heating  and  fiiction  corrections  alone  yielded  M^  =  1.34.  The  last  run  with 
no  corrections  showed  a  Mach  number  of  1.36.  Therefore,  the  mass  loss 
equations  have  a  large  affect  relative  to  the  heating  and  fiiction  corrections. 
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CHAPTER  VI 


CONCLUSION 


BENEFIT 

The  difficulties  in  operating  shock  tubes  in  an  efficient  manner  have 
been  made  evident  in  this  thesis.  Many  problems  arise  in  correctly 
predicting  shock  tube  performance  in  order  to  utilize  its  testing  versatility. 
These  problems  include  burst  diaphragm  interactions,  materials  limits,  and 
shock  wave  attenuation,  to  name  a  few.  These  problems,  however,  are  only 
general  in  nature  and  involve  many  sub-problems  which  contribute  to  the 
overall  understanding  of  shock  tube  characteristics. 

In  order  to  solve  the  attenuation  problem,  several  areas  must  be 
addressed.  These  areas  include,  but  are  not  limited  to:  shock  wave  wall 
friction,  heat  transfer,  real  gas  effects,  and  the  subject  of  this  thesis,  the 
viscous  wall  boundary  layer  interaction  produced  by  passage  of  the  shock 
wave.  This  boundary  layer  interaction  introduces  the  mass  loss  phenomenon 
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to  both  turbulent  and  laminar  flow.  In  turn,  the  mass  loss  affects  the  shock 
wave,  the  contact  surface,  and  the  expansion  wave. 

This  thesis  has  focused  on  a  specific  portion  of  the  larger  problem;  the 
effect  of  the  mass  loss  on  the  shock  wave  for  turbulent  flow  in  a  shock  tube. 
This  work  has  led  to  some  promising  results  throughout  the  range  of  Mach 
numbers  and  including  different  shock  tube  configurations.  The  computer 
code  produced  in  the  process  of  this  study  can  be  applied  readily  to  any 
shock  tube.  The  result  will  improve  shock  tube  characteristic  predictions  and 
allow  for  more  efficient  testing  and  a  better  understanding  of  shock  tube 
performance. 

RECOMMENDATIONS 

The  major  recommendations  consist  of  those  applications  not 
addressed  completely  in  this  study.  The  mass  loss  equations  were  applied 
only  to  the  boundary  layer  between  the  shock  wave  and  the  contact  surface 
because  of  the  fixed  shock,  steady  state  equations.  An  extension  would 
involve  analyzing  the  mass  loss  contribution  by  the  boundary  layer  in  the 
expansion  region  and  the  interface  region. 
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Friction  and  heat  corrections  were  estimated  here  to  provide 
approximate  results.  A  detailed  analysis  of  the  friction  between  the  moving 
shock  wave  and  the  shock  tube  wall  could  be  initiated.  Because  of  the  high 
temperatures  produced  by  shock  tubes,  heating  corrections  between  the 
internal  gas  and  the  shock  tube  wall  should  be  examined.  The  next  step 
would  involve  a  complete  analysis  of  the  real  gas  effects  in  a  shock  tube. 
Depending  on  the  Mach  number,  real  gas  effects  could  account  for  a 
secondary  portion  of  the  prediction  error. 

Lastly,  the  laminar  boundary  layer  could  be  examined.  Mirels 
developed  a  separate  set  of  equations  to  solve  this  regime.  For  this  thesis, 
these  equations  were  dismissed  from  an  analysis  standpoint,  due  to  time 
constraints.  Laminar  flows  may  appear  in  a  large  quantity  of  shock  tubes, 
since  shock  tubes  can  operate  in  the  lower  Mach  number  regions.  Also,  a 
transition  analysis  may  prove  important,  since  a  large  difficulty  exists  in 
actually  determining  when  the  boundary  layer  transition  occurs  inside  the 
shock  tube. 

These  recommendations  are  just  a  few  suggestions  where 
improvements  can  be  made.  Additional  work  can  be  continued  where  this 
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thesis  left  off.  The  results  found  here  have  been  impressive.  The  author 
feels  that  a  significant  step  has  been  made  toward  the  goal  of  better 
efficiency  for  predictions  in  shock  tube  testing. 
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